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Abstract. We present results on number counts of Extremely Red Objects (EROs) in a 2850 arcmin'^ near-infrared 
survey performed in European Large Area ISO Survey (ELAIS) fields a,t K < 17.5. Counts of EROs are extended 
to brighter levels than available previously, giving 0.002 ± O.OOlarcmin"^ at K < 16.5 and consistent numbers 
with literature values at fainter magnitudes. Photometric redshifts from HYPERZ as well as GRASIL model 
SEDs of galaxies imply that our EROs are located in the range z — 0.7 — 1.5, with the bulk of the population 
at 2 ~ 1. Taking advantage of the ISO data in the fields, we use mid-IR detections to constrain the number of 
dusty EROs, and also discuss the superior capabilities of Spitzer Space Telescope to detect dusty EROs. Both the 
mid-IR data and the use of colour-colour diagrammes indicate that at most 10-20% of the EROs in this bright 
regime are dusty starbursting systems. The space density of our EROs, interpreted to be counterparts of local 
> 2 — 31/* massive galaxies at around z ~ 1, is estimated to be « 2 x 10~^ Mpc~^, which is consistent with local 
values. Furthermore, the cumulative number counts at our bright magnitudes are remarkably well fitted by pure 
luminosity evolution models. 

Key words. Galaxies: evolution - Galaxies: formation - Infrared: galaxies - Cosmology: observations - Galaxies: 
starburst - Galaxies: elliptical and lenticular 



1. Introduction 

Extremely Red Objects (EROs, selected for example by 
R — K>5,I — K>A colours) have received much focus 
recently by virtue of their potential as a powerful window 
into the galaxy formation era. By and large, the majority 
of EROs have been incorporated into a bimodal popula- 
tion, where the extreme red colours are attributed either 
to old passively evolving distant {z > 1) elliptical galaxies 
or to extremely dust reddened starburst galaxies (see e.g. 
recent papers by Cimatti et al. 120031 Wold et al. 120031 Yan 
& Thompson '^003', Takata et al. "211031, Daddi et al. 00021 
Small et al. 2002, Roche et al. 2002, and references therein 
for earlier pioneering ERO surveys). It is the class of aged 
early type galaxies whose number densities provide the 
strongest constraints on models of galaxy evolution. On 
the other hand, the dusty EROs could be related to the 
(ultra) luminous IR-galaxies producing the bulk of the to- 
tal energy in the Universe since the recombination era (see 
e.g. Elbaz & Cesarskv 120031 for a review). 



The dusty ERO population contaminates the other- 
wise, in principle, clean samples of high redshift massive 
ellipticals. This is the reason why the separation of ellip- 
tical EROs from dusty EROs, using photometric, spectro- 
scopic, and morphological analyses, has received the main 
attention of ERO studies recently. In general, a large num- 
ber of massive galaxies at redshifts of unity and over would 
support a traditional monolithic collapse scenario where 
galaxies form at high redshift in a single collapse and then 
evolve passively (e.g. Eggen et al. 119621; Larson 11975(1 . 
Significantly smaller numbers of massive high- 2; ellipticals 
on the other hand fit predictions from the hierarchical as- 
sembly scenarios, where massive galaxies are formed only 
relatively recently from the merging of smaller units (e.g. 
White & Rees [W751 White & Frenk [TMTl Somerville & 
Primack 1999: Cole et aLlMKl. 



Send offprint requests to: P. Vaisanen 
e-mail: pvaisaneOeso . org 



Since EROs are extremely faint optically (i? > 24), 
their redshifts and spectral and morphological proper- 
ties have remained largely unknown until very recently. 
Cimatti et al. (|20D2i) carried out a VLT spectroscopic sur- 
vey of 45 EROs with R-K > 5 and A% < 19.2. They iden- 
tified approximately 1/3 of the EROs as old elliptical sys- 
tems, 1/3 as dusty starburst galaxies, while 1/3 remained 
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unidentified. The mean redshift was found to be 2 ~ 1.0 
for both populations. Yan, Thompson, & Soifer H2004() on 
the other hand find a small 10-15% fraction of isolated 
passive systems. Yan & Thompson (2003), Cimatti et al. 
iPTHSjl . Moustakas et al. (2004), and Gilbank et al. (:200g| ) 
have used HST morphologies to differentiate between the 
classes: the results have further complicated the picture, 
as a large fraction (25-65%) of EROs seem to be disks 
(though see also Moriondo et al. I2()00|l . The expectation 
originally had been a clearer distinction into spheroids 
on one side, and irregular and interacting types on the 
other, expected for extremely dusty star forming galaxies. 
The relative fractions of different types of EROs thus still 
remain uncertain (due to for example different selection 
criteria used) and the results of testing surface densities 
against galaxy formation scenarios are inconclusive. 

In addition to spectroscopic and morphological meth- 
ods mentioned above, one may separate the old elliptical 
and dusty EROs by means of longer wavelength observa- 
tions: any ERO detected at mid-IR to radio wavelengths 
should belong to the dusty population. Systematically this 
has been attempted by Mohan et al. (|2002|) in the sub- 
mm to radio and Small et al. H2002|l in the radio. Similar 
fractions (albeit with wide spread) as referred to above of 
dusty EROs were found in the latter study, while Mohan 
et al. (.2002.1 find a much lower fraction. 



In this paper we attempt to separate EROs, for the 
first time, based on their mid-IR properties. As with far- 
IR to radio methods, this is a clear-cut definition, since the 
difference of ellipticals and starbursting galaxies is very 
large - distant evolved ellipticals would not be detected 
with ISO whereas dusty EROs should have strong mid-IR 
flux. Also, mid-IR avoids the identification problems of e.g. 
large sub- mm beams. We have surveyed an ISO/ELAIS 
(European Large Area ISO Survey; Oliver et al. 12000)) field 
in the near-IR and matched the results with an optical 
dataset. We make use of the newly available Final ELAIS 
Catalogue (Rowan- Robinson et al. 12004)) . In this paper 
we consider the surface densities of EROs resulting from 
a wide-field relatively shallow J and if-band survey. We 
shall present the results from a deeper NIR survey around 
faint mid-IR detections in subsequent papers (Vaisanen & 
Johansson 120041 Johansson, Vaisanen & Vaccari 2004, in 
prep.). 

The structure of this paper is as follows. The observa- 
tional data are presented in §2, and in §3 we define and 
extract the ERO sample from the photometric catalogues. 
In §4 we discuss the dusty vs. evolved ERO separation and 
in §5 the SEDs of the detected EROs. The resulting sur- 
face densities of EROs are presented in §6. We assume 
throughout this paper a flat (fig — 1) cosmology with 
VLm = 0.3, = 0.7 and = 70 kms^^Mpc-i. 



2. Observations 

2.1. Near-Infrared data 

The observations presented here were conducted in 
the ELAIS fields Nl and N2, centered at {a,S) = 
16h09m00s, 54°40'00" and (a, S) = 16h36m00s, 41°06'00" 
J2000.0 respectively. 

The near-IR survey performed with the Mt. Hopkins 
1.2-m telescope, reaching limiting magnitudes of approx- 
imately J = 19.3 and K = 17.5, is fully presented in 
Vaisanen et al. (j2000i hereafter VOO). The data cover a 
total of one square degree, although a slightly smaller area 
is used here (2850 arcmin^) after some edges and higher 
noise areas were excluded and, more importantly, only ar- 
eas with full coverage in both bands were considered. The 
detector used was a 256 x 256 InSb array with 1.2" pixels. 

2.2. Optical data 

We obtained optical photometric data from the Isaac 
Newton Telescope (INT) Wide Field Survey (WFS). 
The WFS data are publicly available on the Cambridge 
Astronomical Survey Unit (CASU) homepage^ For a re- 
view on the INT Wide Field Survey Project and instru- 
ment characteristics, see McMahon et al. ()2()()l)l . The data 
consist of U, g' ,r' ,i' , Z band photometry (Vega-based) in 
the Nl region and g',r',i',Z photometry in the N2 re- 
gion (the apostrophes are dropped henceforward for clar- 
ity). The WFS bandpasses are similar to the SDSS (Sloan 
Digital Sky Survey) filters (Fukugita et al. ITM5|) . The 
WFS webpage^ gives colour transformations between the 
WFS filters and those of the standard Johnson-Cousins 
system (Landolt EMS) : r - i? = 0.275(i? -/)+ 0.008 and 
i — I = -|-0.211(i? — J) for the bands most used in this 
paper. The nominal Scr detection limits for a 1" seeing 
are g « 25.0, r « 24.1, i « 23.2,and Z « 22.0 (the WFS 
page). For further details on the pipeline processing of 
INT wide field survey data consult Irwin & Lewis (,20(TT]) . 
and Gonzalez-Solares et al. ()2004)) for more in-depth dis- 
cussion of the WFS data in particular in ELAIS fields. We 
also obtained raw i?-band images of the ELAIS Nl and N2 
regions from the archives. After reduction we used them 
for object identification - only the pipeline processed cat- 
alogues were used for photometry however. 

2.3. Astrometry and photometry 

The INT astrometry was derived using Guide Star Catalog 
(GSC) stars, which results in an external astrometric ac- 
curacy of 0.5" — 1.0" (Irwin & Lewis i2001) . Astrometric ac- 
curacy in the Mt. Hopkins data ranges between 0.5 — 1.5", 
and was calibrated using GSC and US Naval Observatory 
(USNO) catalogues. 

All our near-IR photometry is performed using the 
SExtractor software (v. 2. 2.1 and v. 2. 3; Berlin & Arnouts 

^ http://archive.ast.cam.ac.uk/ 

^ http:/ /www.ast. cam.ac.uk/~wfcsur/index.php 
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I1996|l . Our NIR photometry is explained in VOO. As dis- 
cussed therein, we used total magnitudes as given by the 
SExtractor 'BEST'-magnitude. However, since the optical 
WFS photometry was given only in 2.4" diameter aper- 
ture magnitudes, we re-did all the old photometry with 
matching apertures. It is crucial that the different obser- 
vations map the same region of the source when construct- 
ing colour indices; near-IR magnitudes not associated with 
colours are given as total magnitudes in this paper. 

The archive optical INT WFS data in the ELAIS re- 
gions come as fully calibrated photometric source cata- 
logues (Irwin & Lewis 200T Gonzalez- Solares et al.|2004). 
No further processing is done for the catalogue except 
merging of multiple and/or nearby sources, as described 
below in Section 

3. Construction of the ERO samples 

3.1. Definition of EROs 

Numerous different selection criteria have been defined for 
EROs, including R-K>6, R-K> 5.3, R - K > 5, 
I ~ K > A with iiT-magnitude upper limits from 18 to 
21 mag. All these criteria are designed for selecting early 
type galaxies at z > 1. In this paper we use the following 
definition for EROs: r — K > 5.5 and/or i — K > 4.4. Our 
limits result from colour transformations given above in 
Section for typical colours ofi? — / »1.4 — 2.0of our 
sources, and from the desire to have limits corresponding 
to the commonly used R — K>5, I~K>4: selections 
to aid comparisons to other surveys. In addition, we will 
check any results with r — K > 5.8, i.e. R — K > 5.3, which 
is also often used. 

Furthermore, we calculated r — R and i — I colours us- 
ing SEDs from the GRASIL hbrary (Silva et al.Ll998) and 
verified the r — Rk, 0.5 and i — I k, 0.4 colours to match 
well those of ellipticals at 2 ^ 1. Naturally an exact com- 
parison or transformation between different galaxy colour 
selection criteria is redshift and galaxy SED dependent, 
and we note that since the model colours that different 
authors use vary, there might easily be 0.1-0.3 mag differ- 
ences in the colours at z 1. Elliptical galaxies become 
EROs atzsal.l — 1.2 with the criteria and models we 
used. 

Fig-fflshows r — K and i ~ K model colours of several 
representative galaxies against redshift with the ERO cri- 
teria included. Model SEDs are adopted from the GRASIL 
library (Silva et al. [1998; model SEDs and the GRASIL 
code are available online'^). Ordinary spirals (dotted line) 
never reach the red colours of EROs, while both ellipticals 
(solid) and reddened starbursts (dash-dot) become EROs 
when seen beyond z ~ 1. For comparison, the colour of the 
prototype dusty ERO HRIO is also plotted (dashed curve) 
as a function of redshift. The colour of the HRIO model 
is due to extreme dust extinction fSilva lTMSp . The lowest 
panel shows the flux ratio fi5fj,m/ f2.2fj.m - the degeneracy 

^ http:/ /web. pd.astro.it/granato/ 
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Fig. 1. Observed r — K and i — K and MIR/NIR colours 
of representative galaxies are calculated from model SEDs 
of GRASIL code (Silva et al. [TMH| . An evolving Sb spiral 
(dotted line) is plotted along with an evolving elliptical 
(solid curve). Non-evolving models fitted to the observed 
rest-frame SEDs of the prototypical starburst M82 (dash- 
dot) and HRIO ERO (dashed), are also plotted. Solid hor- 
izontal lines show the colour definition of EROs in this 
paper. The dash-triple-dot curve in the lowest panel cor- 
responds to the ULIRG ARP220, which gives the most 
extreme MIR/NIR colours. 



in the red colours of old ellipticals and dusty starbursts is 
clearly broken, as will be discussed later in more detail. 

3.2. Matching of optical and NIR data 

The optical counterparts of near-IR sources in the J and 
if-band matched Mt. Hopkins catalogue were extracted 
from the INT WFS catalogue using 3" search radius. More 
than 80% of the matches are found within « 1.5" separa- 
tions, which is reasonable given the astrometric accuracies 
and pixel sizes. Approximately 4% of the 6600 NIR sources 
in the Mt. Hopkins catalogue were not matched with an op- 
tical source. We checked these individually: the great ma- 
jority are close to WFS CCD edges or bright stars, or were 
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some remaining NIR frame cosmic rays, etc., and though a 
handful of these may be genuine very red objects (though 
not necessarily EROs) we decided to conservatively ex- 
clude all of these from the final source list. 

Since the WFS catalogue contains multiple detections 
of the same source (due to overlapping CCD frames), opti- 
cal counterparts within 0.5" of each other were averaged. 
After this purging, the remaining multiple matches within 
1.5" were summed up since the corresponding NIR cata- 
logue would not have resolved them as individual sources, 
and then the brightest source was selected if there still 
were multiple counterparts available. Approximately 18% 
of the original WFS optical catalogue closest matches were 
affected by the purging, though less than 4% by more than 
0.1 magnitudes. 

3.3. Star vs. galaxy separation 

We then proceeded to separating stars from galaxies in the 
surveys. Stellarity indices were available from both the 
NIR catalogue (the SExtractor CLASS parameter) and 
the WFS catalogue, where a flag defines galaxies, definite 
stars, and various degrees of uncertain stellarities. The 
WFS classification is in principle more useful for us here, 
since it goes deeper than the NIR data. Note that the 
star-galaxy separation is different from that conducted in 
VOO since we now have a much deeper optical catalogue 
available. 

However, for the EROs we are interested in, the red- 
dest and faintest objects in the catalogue, we find that a 
colour separation works best. This can be seen in Fig.|21 
which shows the full catalogue m r — i vs. r — K with 
the morphological classification indicated: stars are over- 
plotted as small crosses in the left panel. A separating 
line, adjusted experimentally to maximally distinguish the 
main concentrations of stellar and extended sources in the 
figure, oir ~ K — 2.16(r — i) + 1.35 is drawn. This colour- 
colour diagramme separates stars very well from galaxies; 
it is important to note that there is virtually no overlap 
between the stellar sources and galaxies in the region of 
interest aX r — K > 5. The r — i colour, in fact, is very 
important in the separation of red stars. Very low mass 
and cool stars (of the L spectral type in particular) have 
R — K and J — K colours closely mimicking those of extra- 
galactic EROs - for these stars, however, R — I (or r — i) 
always stays above « 2 (see e.g. Chabrier et al. 120001 Cruz 
et al. OnnSji . On the other hand, we find no galaxy models 
resulting in r — i > 2 colours (GRASIL code used; see also 
eg. Fugukita et al. I1995|l . All extremely red objects with 
r — i > 2 can safely be discarded as stars. 

Representative galaxy colours calculated from 
GRASIL models (Silva et al. I1998|l are overplotted in 
the right panel. In fact, the elliptical model does overlap 
with the r — j ~ 3 extremely red stars, but only at z > 4. 
This is a potential concern with deeper FRO surveys, but 
at our brighter magnitudes such distant sources are not 
expected to be seen. 



We thus use a combination of methods. We define as 
stellar all objects having stellar colours according to the 
above limit and those not having a galaxy morphology set 
by the WFS survey. The end result of this is that brighter 
sources are preferentially separated by morphology and 
fainter ones (especially red sources) by colour. The star 
counts were already compared to the SKY model predic- 
tions of Cohen 1)1994(1 in VOO, and found to fit well model 
predictions for the corresponding fields. In the regime 
K < 17.5, 45% oi r — K > 5.5 sources were classified 
as stars. In summary, it should be stressed that in the 
case of bright surveys of EROs, the contamination from 
red late type (typically L-type) stars is considerable when 
using the R — K based selections of EROs. 

3.4. Extraction of EROs 

We searched for EROs from the Mt. Hopkins galaxy cat- 
alogue according to the colour definitions given above in 
Section im Only NIR detections at 5a level and over were 
considered. We then went through the resulting list check- 
ing our NIR maps as well as the WFS CCD images for 
any obvious spurious objects (for example, anomalously 
low optical magnitudes were found in cases when sources 
fell close to gaps in the CCD frames, or near bright stars 
- these were excluded) . Ultimately there were 50 EROs in 
the matched Mt. Hopkins catalogue using the r ~ K crite- 
rion and 21 using the i — K limit - 17 EROs are common 
to lists resulting from both selection criteria. This makes 
54 EROs in total. All the EROs have an i-band detection, 
while there were 4 EROs which have only upper limits 
in the r-band. It is obvious that the r and i band based 
selection criteria for EROs are not equivalent. There are 
more than twice as much r — K > 5.5 selected EROs than 
i — K > 4.4: selected ones. 

The photometry for all EROs is given in Tableland 
the resulting total numbers of verified EROs are summa- 
rized in Table El along with some other survey character- 
istics. The colour-magnitude plot is shown in Fig. |31 

We wish to point out, that in order to be very con- 
servative in limiting the number of spurious detections in 
final ERO lists, as well as to be able to take advantage of 
the J — K colour in classifiying EROs, we required both 
a J and a K detection for all considered objects from the 
survey. While the 5cr detection requirement was applied 
only to the /^-band, this nevertheless might exclude some 
genuine EROs which are very red in J — K colour. At our 
faintest levels of K > 17.5, we do not expect to detect any 
J — K > 2.5 EROs, and already dX K >n the complete- 
ness oi J ~ K > 2 EROs would not be quite as high as 
J -K <2 EROs. 

4. Separation of dusty EROs 

Since the Extremely Red Objects can be divided into two 
broad classes, the populations should be differentiated be- 
fore surface density comparisons to detailed galaxy for- 
mation models are attempted. A prime motivation of this 
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Fig. 2. Observed r — i vs. r — K colours of our survey. Stellar sources, as defined by the WFS classification are 
overplotted as small crosses in the left panel. It is seen that this colour-colour diagramme does a very good job in 
separating stars from galaxies: in particular there is virtually no overlap in the region of interest, at r — K > 5. Any 
"ERO" with r — i > 2 is classified as a star. The right panel shows all the galaxy EROs of the Mt. Hopkins survey 
selected with r — K > 5.5 and several GRASIL model SEDs overplotted: the lowest solid curve is the SED of an Sb 
galaxy, the dashed curve reaching ERO regime is that of a starburst, and the highest curve an elliptical. The curves 
are plotted up to redshift of z = 2 (see text and Sect. 



study is to attempt the separation using the mid-IR data. 
We shall also perform the separation using other estab- 
lished methods, and discuss the differences. 

4.1. Using mid-IR data 

Since only dusty EROs are expected to show up in the 
mid-IR ELAIS survey, the question to answer is, what 
fraction of the EROs are detected in the mid-IR to a given 
flux limitl 

Among the 54 Mt. Hopkins EROs we find only one 
match with a 15 /xm ISOCAM source. In the ELAIS band- 
merged catalogue (Rowan- Robinson et al. l2004|l the source 
is a candidate hyper luminous IR-galaxy based on a pho- 
tometric redshift of z ~ 1.0, which agrees with our inde- 
pendent photometric redshift determination oi z — 0.9. 
(We note that using the same near-IR data-set but with 
only if -band considered, another ISOCAM detected ERO 
is found in Rowan- Robinson et al. H2004|l ). Thus, at face 
value, the fraction of dusty EROs seems to be insignicant. 

However, it is crucial to take into account detection 
limits. The main hindrance to the full power of the mid- 
IR separation with our data-set is the relative shallow- 
ness of the ELAIS survey. All the Mt. Hopkins EROs have 
if = 16- 18 mag. Comparing to GRASIL models (Fig. HI), 
a starburst galaxy with M82-type SED and K — \7 ap- 
parent magnitude should just have been detected over the 
0.7 mJy ISOCAM flux limit. In comparison, to detect 



the Rayleigh-Jeans tail of ellipticals would have required 
approximately if = 13 brightness. We note that M82- 
type galaxy with K — 17 mag dX z — 0.7 and z = 1.0 
translates to star formation rates SFR « 100 and 200 
M0yr~^, respectively, as calculated from total IR lumi- 
nosity iiR(3 — 1000/J,m) with relation given e.g. in Mann 
et al. (20ll|) and using the GRASIL model SED. The cor- 
responding IR luminosities are log^Lig) > 11.7. 

There are 23 if < 17 EROs in the sample, out of which 
only one (4%) is detected in the mid-IR. As shown above, 
the ELAIS data is only sensitive to luminous IR galaxies 
at the expected redshifts of EROs (z > 0.7). Thus, we can 
state that very strong starbursts [SFR > 100 MQyr~^ 
make up only a small fraction, less than 10%, of counter- 
parts to bright EROs. For detection of more modest dusty 
galaxies other methods have to be used. 

It is relevant to note that the in-orbit Spitzer mission 
will definitely find large numbers of dusty EROs, including 
more modest ones than above: Using the M82 SED once 
again as an example, sampling dusty if « 20 EROs at 
z « 1 would mean probing dusty starforming galaxies of 
SFR lOMQyr"^. The expected flux densities of such 
objects with IRAC 8/ini and MIPS 24/ini bands would 
be of the order of 10 and 150 /iJy, respectively, which 
are easily reached in a few minutes of integration time. 
For example, assuming 10 minute per pixel integrations, 
a 10 hour survey to the mentioned 5a depths would cover 
~ 800 arcmin^ in both bands. This means detecting 100- 
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Fig. 3. The iiT-band magnitude vs. the r ~ K and i ~ K colours, in left and right panels, respectively. All objects 
coming through the catalogue construction as galaxies are plotted as open circles. Individually verified EROs at > So- 
level are overlaid with solid circles. Since the survey was done in ELAIS fields, we also overplot all matched ISOCAM 
detections as solid triangles. Note that the total iC-magnitude is used on the x-axis, whereas the colour is calculated 
with matching small apertures. 



200 MIR counterparts of dusty EROs, assuming average 
ERO surface densities and an ad hoc 50% fraction of dusty 
EROs. 



4.2. Colour-colour selections 

Figure 0] shows the colour-colour separation scheme of 
EROs by Pozzetti & Mannucci (|2000)1 adopted for our r, i 
filters. The idea is that the J — K separates the EROs at 
z > 1 to bluer early type galaxies (where the large optical 
to near-IR colour is due to the 4000A break), and to the 
redder (in J — K) dusty EROs which have a a smoother 
SED from optical to near-IR. In the r — K plot we find 
9/50 r — K selected EROs in the dusty starburst side of 
the indicator, and similarly 3/16 of the i — K EROs. In 
the i — K plot many r — K EROs fall under the i — K — 4.4 
line, and the dusty percentages are higher for the remain- 
ing EROs, 30 - 50%. The one mid-IR ERO is overplotted 
in the figure, and it in fact falls on the elliptical side of 
the division. However, given the typical photometric er- 
rors of the EROs (the error bars on the mid-IR ERO are 
representative) it is seen that at least 1/3 of the EROs fall 
statistically on the dividing line. Moreover, while the ellip- 
tical vs. starburst separating line is equivalent to Pozzetti 
& Mannucci, it is only defined to work for R — K > 5.3. If 
we thus select EROs at r — K > 5.8, nearly all galaxies on 
the starburst side of the left panel in Fig. 0] fall out, and 
the elliptical fraction becomes close to 90%. Finally, our 
Mt. Hopkins sample is 1-2 magnitudes brighter than other 



typical ERO samples, and it is therefore not certain how 
accurate the separation should be. 

We also checked the colour-colour separation scheme 
presented by Bergstrom & Wiklind (|2004|l using R — J 
vs. J — K colours (see their Fig. 8). The result is 85% 
ellipticals regardless of whether r — Koii — K ERO crite- 
rion (or both) is used, totally consistent with the fractions 
emerging from the Pozzetti & Mannucci R — K vs. J — K 
method. 

We thus conclude that using the various J — K vs. 
optical-NIR colour separation schemes approximately 80% 
of our EROs appear ellipticals. Only when using Pozzetti 
& Mannucci i — K based separation for i — K> 4.4 EROs, 
is the percentage lower, « 60%. We might be missing the 
very reddest J — K EROs, i.e. some dusty EROs, because 
of the selection of the sample in both bands. At K < 17, 
where there should be no bias in J— if colour, the fractions 
of early type EROs range between 65-80%, depending on 
the optical band used. 

5. Photometric redshifts and SEDs 

We calculated photometric redshifts for the EROs using 
the HYPERZ code pWijl which fits GISSEL98 synthetic 
spectra (Bruzual & Chariot 119931 BC hereafter) galaxy 
templates to photometric data points. Figure shows 
the photometric redshift distributions of our sample. The 
mean photometric redshift is z = 0.94 ± 0.49, excluding 
the outliers at z ^ 4. However, it is clear that subgroups 
are involved, though the gap at z = 1 — 1.2 is difficult to 
understand. We also note that approximately half of the 
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fits result in a confidence better than > 85%, so the results 
should be taken with some caution. On the other hand, 
the overall distribution, average redshifts, or properties of 
sub-groups discussed below do not change if we consider 
only those EROs with fits of > 85% confidence. 

The 21 i — K selected EROs, the white regions in the 
histogram, are more evenly distributed and at higher z 
(average z w 1.4) than r~K EROs. The 17 objects which 
are EROs with both selection criteria cover quite evenly 
the range z = 0.6 — 1.8 (average at z « 1.3). The EROs 
#51-54 in TableUare all between z = 1.4 - 1.9. 

Most notably, however, those EROs that are not EROs 
with the i — K > 4.4 criteria, show a significantly nar- 
rower redshift distribution (hatched region of histrogram) 
with z = 0.73 ± 0.06. Moreover, these objects constitute 
more than half, 33/54, of our total ERO sample - we re- 
turn to these below in Section 16.21 Taken together, the 
photometric redshifts thus suggest that r — K > 5.5 se- 
lects more nearby systems at redshift unity and below, 
and i — K > 4.4 includes a wider range of EROs. 

Virtually all, 50/54, of the best-fit HYPERZ BC SEDs 
are either starburst spectra of age 0.4 to 2 Gyr, or ellipti- 
cal spectra at several Gyr - these SEDs are nearly iden- 
tical (see Fig. EJ. In fact, in a recent paper by Pierini et 
al. (|2Q04) it was proposed that post-starburst ellipticals, 
forming in a short bursts during the period where EROs 
are observed (between 1 < z < 2), are a neglected con- 
stituent of ERO populations. It should be stressed that ba- 
sically all of the best-fit HYPERZ 'starbursts' are evolved 
results of instantenious bursts and not dusty starforming 
galaxies. 



20 



15 



10 



Photometric redshifts 
of Mt. Hopkins EROs 




12 3 4 

redshift 

Fig. 5. Photometric redshifts as calculated by HYPERZ. 
All 54 Mt. Hopkins EROs are plotted, whatever the confi- 
dence of the fit (half have > 85%). r — K > 5.5 selected 
EROs which are not EROs with i — K > 4.4 selection are 
hatched, the empty regions of the histogram thus show 
thei-K > 4.4 EROs. 
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Fig. 6. Averaged and if-band normalized SEDs of i — 
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GRASIL M82 model and elliptical and evolved starburst 
SEDs from BC are also plotted at indicated redshifts. 

6. Discussion 

6.1. ERO counts and surface densities 

The numbers of detected EROs are given in Tabled and 
the corresponding cumulative number counts are plotted 
in Figd The numbers plotted are the total ERO counts: 
based on the previous discussion, the fraction of dusty 
EROs is low (~ 15%) and does not have significant effect 
even if subtracted from the counts. 

Our wide survey from Mt. Hopkins extends the ERO 
number counts to brighter magnitudes than observed be- 
fore. We report the first ERO counts a,t K < 16.5, where 
the r — K > 5.5 (ie. R — K > 5) selection yields a surface 
density of 0.002 ± O.OOlarcmin^, and the i - K > 4.4 (ie. 
/ - X > 4) selection w 50% of this value. At ii' < 17 mag- 
nitudes we arrive at 0.007 ± 0.002arcmin^ oir — K EROs. 
Beyond K > 17 mag the Mt. Hopkins survey is starting to 
be affected by incompleteness, and based on simulations 
performed in VOO for the same dataset, we estimate a fac- 
tor of 1.5 correction to the K < 17.5 cumulative ERO 
density. The resulting r — K selected surface density at 
K < 17.5 is 0.018 ± O.OOSarcmin^. 

We compare these numbers to the only two other 
field ERO surveys with enough sky coverage to reach 
these bright magnitude levels: The largest ERO survey 
to date Daddi et al. H2()()()al with 701 arcmin^), yielded 
0.003 ± 0.002arcmin2 at < 17 and 0.02 ± O.OOGarcmin^ 
at K < 17.5. These are totally consistent with our result. 

The Yan & Thompson ^^MM 409 arcmin^) HST ERO 
counts were selected using I — K > 4.0. At X < 17 they 
agree with our corresponding i — K > 4.4 selection within 



the errors. At i^T < 17.5 there is a clear discrepancy, how- 
ever. It is our r — K selected EROs which are in closer 
agreement with the Yan & Thompson ERO counts - we 
systematically find a factor of ~ 2 less i — K selected 
EROs than r — K selected ones. Especially at the faintest 
bin there appears to either be more incompleteness than 
we expect, or the i — K > 4.4 selection does not pick up 
as many bright K ^ \7 galaxies at z ^ 1 than the r — K 
selection does. Indeed, the discussion on "excess" r — K 
EROs in the next Section 16.21 and the photometric red- 
shifts derived for the EROs hint at the latter possibility. 

6.2. Evolved ellipticals at z ~ 1 

The strength, and motivation, of ERO searches has been 
to look for the most massive galaxies at redshifts of unity 
and over (see e.g. Saracco et al. l20U3fl . since it is exactly 
these which place tight constraints on galaxy formation 
scenarios. The dusty vs. early type ERO discriminators 
showed that most of our ERO sample is consistent with 
being early type galaxies. However, we argue that the 
clearest population of ellipticals comes from the class of 
r — K > 5.5 EROs which have r — K < 4.4 - these consti- 
tute 60% of all our EROs. Fig. IHlshows an i — if vs. r — K 
plot, where these "excess" EROs populate the upper right 
quadrant (see also Figs. 01 and inj. 

Intuitively, the upper left quadrant galaxies should be 
sources at a redshift where a significant 4000 A break falls 
right in between the r and i bands: i.e. ellipticals or very 
early type spirals at a redshift z ~ 0.8. We calculated 
numerous GRASIL models, and it turns out to be fairly 
difficult to obtain colours in this ERO regime. Any signifi- 
cant on-going star- formation drops the r — K colour down, 
as noted also by e.g. Yan & Thompson 1)20031 see their 
Fig. 9.). The effect is seen in the solid curve of the right 
panel of Fig. |S1 the elliptical model has a formation red- 
shift of = 3 and at redshifts of z > 1.5 the effects of the 
starburst turn the r — K colour sharply down. Significant 
amounts of extinction, including edge-on spirals, do not 
bring models to the r — K> 5.5, i — K < 4.4 region either: 
the dotted line depicts the extremely dusty HRIO model 
as an example. 

One gets closest to the r — K > 5.5, z — if ~ 4 re- 
gion by merely redshifting an old present-day elliptical to 
Sa galaxy to the appropriate redshift. This is shown as 
the dashed curve, a 13 Gyr old Sa spiral; any > 10 Gyr 
old elliptical produces a similar result. While this is an 
unphysical model, it serves to point out that while pure 
reddening is unlikely to reach this part of the colour-colour 
diagramme, an old stellar population does that. We thus 
conclude that EROs which are selected by r — if > 5.5 
(i.e. R - K > 5) which are not EROs by z - if > 4.4 
(i.e. i — if > 4) arc mostly early type galaxies. In fact, 
at z — if < 4.2 there should not be any contamination 
from dusty galaxies (compare the vertical dotted line to 
the dusty model curves). 
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We also plot a sample of averaged SEDs of our EROs 
in Fig. El The "excess" EROs wfth i - K < 4.2 are plot- 
ted with squares and are well fit by very evolved galaxy 
SEDs, as just discussed. The critical observed band is the 
i band, where the difference between elliptical and dusty 
starbursts (M82 plotted) is the largest at the redshifts of 
our targets. Note that the (/-band is biased to the bluest 
sources, since the average is calculated from those EROs 
detected in the corresponding band, and only 25% of the 
sources have a g-band detection. 

The r — K and i — K selections alone are quite differ- 
ent (see e.g. Scodeggio & Silva l2UUn|l . At least with bright 
ERO surveys such as ours, it seems that r — K selects more 
passive systems. This is consistent with spectroscopic fol- 
lowup of R ~ K > 5.3 surveys finding significant amounts 
of ellipticals (e.g. Cimatti et al. .2002^1 and those following 
wp I — K > 4 selections finding much more star formation 
(e.g. Yan et al.|lIH|). 

6.2.1. Number densities 

There are 33 "excess" r-K EROs with i-K < 4.4. They 
appear to be at redshifts z = 0.7 — 1.2: BC SED fitting 
favours redshifts below unity (Fig. O , while the GRASIL 
models suggest redshifts at z « 1 — 1.2 (eg. Fig.lHJ. These 
EROs have typical magnitudes of K ^ 17.2, which would 



make them at least ~ 2L* galaxies, taking into account 
passive evolution since that redshift. Very simply, assum- 
ing that this group of EROs consists of ellipticals and is 
constrained in the mentioned redshift range we derive a 
co-moving volume of 2 x 10^ Mpc^ for the Mt. Hopkins 
survey area, and thus a space density of ~ 2 x 10^^ 
Mpc~^. Taking i^-band LF parameters from Kochanek et 
al. H2001|l . this is a significant part of the expected density 
of w 8 X 10~^ Mpc~^ for local > 2L* early type galax- 
ies. Given the incompleteness in our survey a,t K > 17, 
and that we did not include ellipticals a,t i — K > 4.4 in 
this calculation these ERO counts are a very conservative 
lower limit, and thus indicate that a major part of massive 
present day ^ 2 — 3L* ellipticals were in place at 2; ^ 1. 

The space density just derived is intriguingly close to 
those estimated for the higher redshift sub-mm popula- 
tion, which are suspected to be the most luminous star- 
bursts, perhaps the progenitors of present day massive el- 
lipticals (see e.g. discussion in Scott et al. I2002|l . 

6.3. Constraints on galaxy formation scenarios 

Figure |7| plots our ERO counts along with several other 
recent surveys and pure luminosity evolution (PLE) mod- 
els of Daddi et al. H2000a|l . All curves depict ellipticals 
with a r = 0.1 Gyr initial star burst, and passive evolution 
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Fig. 8. i — K vs. r — K colour-colour magnitude diagramme of the verified EROs is plotted in the left panel. The 
mid-IR detected EROs is overplotted as a large solid symbol. Several GRASIL models are plotted in both panels: the 
models are drawn until redshift z = 2, and the z = 1 locations are indicated by small solid circles along the curves. 
The reddening direction is indicated in the left panel by the diagonal arrow which correcponds to Ay = 1 extinction. 



afterwards. The different lines have formation redshifts of 
Zform = 2.4,3.0, and 10. All the models use a Marzke et 
al. ([195^ LP, except the uppermost dashed line, which 
instead uses a 2MASS based LF (Kochanek et al. [2001 ) 
resulting in a higher normalization by approximately a 
factor of two. 

We stress that all EROs are counted in the figure, i.e. 
no corrections for the dusty population have been made. 
As was seen previously, this is perfectly reasonable for our 
very bright sample of EROs. However, even a very large 
correction (not supported by our observations) of 50 % 
of the EROs being in fact dusty, would not change the 
conclusion significantly. At somewhat fainter magnitudes, 
recent studies have shown (Cimatti et al. 200?, Small et 
al. l2002|l that the percentage of dusty EROs is likely some- 
where around 30-60%. 

The PLE model using a formation redshift of Zform = 3 
fits all the data remarkably well. Pormation redshifts 
slightly lower than this also fit well our bright ERO counts, 
but start underpredicting the numbers at fainter magni- 
tudes. On the other hand, corrections for dusty galaxies 
might be larger at the fainter levels (see eg. Smith et al. 
I2002|l . However, models with Zform < 2.4 start underpre- 
dicting the counts at all levels, especially when using the 
R — K > 5.3 selection criteria. The very highest formation 
redshifts of Zform = 10 and over, predict steeper number 
counts and arc more difficult to fit to both the bright and 
faint end of ERO counts with the same normalization. The 
figure shows two different LPs giving a factor of ^ 2 differ- 
ence. Note also that eventhough the R~ K > 5.3 criterion 
eliminates nearly all the (possibly) lower redshift z ~ 0.8 



EROs, there is no significant difference in the fits of PLE 
models in the right and left panels of Pig. [3 

Using the 2MASS normalized LP it is in fact slightly 
more difficult to fit both R - K > h and R~ K > 5.3 
counts simultaneously. The R~ K > 5 counts are overpre- 
dicted by a factor of 1.5 — 3 with all except the very lowest 
Zform = 2.2—2.4 formation redshifts. Good fits in the range 
K — 16.5 — 18.5 are again acquired with Zform > 4, but the 
faintest counts are overpredictcd. Roche et al. H2002I2003|I 
find their 2MASS-LP based PLE models significantly over- 
predicting I — K > 4 ERO counts, and show that certain 
amount of merging and density evolution fits the counts 
much better in particular in the fainter i^T > 20 regime. 

That a range of formation redshifts are necessary to 
model EROs is seen for example by comparing Pigs. |H1 
and (see also Cimatti et al. I2()()3|l : though models with 
formation redshifts of z « 3 fit well the counts, the colours 
produced by such a scenario are not red enough for many 
EROs found in this survey, and others. Very red colours 
oi r — K > 6.5 or i — K > 5 can not be produced without 
pushing Zform closer to 10 (unless all such extreme EROs 
are heavily reddened starbursts). 

We do not investigate hierarchical formation scenarios 
in more detail here. We merely point out that predictions 
from some recent models (Cole et al. l2UUUI as presented in 
Smith et al. l2002|l . fall short an order of magnitude in the 
numbers of EROs in the range K = 17 — 20. See also e.g. 
discussion by Martini (PTOTIl and Pirth et al. H2002|l . It is 
important to realize that a drastic increase in the fraction 
of dusty EROs does not make the hierarchical models fit 
the ERO counts any better: they predict too few EROs of 
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all kinds. On the other hand large fractions of dusty EROs 
would have to be subtracted from PLE models including 
only ellipticals, making all the Daddi et al. PLE models 
used above to overpredict the counts by factors of 1.5 — 3. 

It is concluded that the PLE models do give remark- 
ably good fits to the brightest ERO counts. Formation 
redshifts around z ~ 3 are favoured - however, by alter- 
ing the details of the models, formation redshifts between 
z = 2 — 10 are also consistent. Moreover, a range in the 
formation era of ellipticals is suggested by the range in 
ERO colours. 

7. Summary 

We have searched for EROs in a near-IR survey performed 
in ELAIS fields, using r - K > 5.5, r - K > 5.8, and 
i — K > AA colour criteria. These are equivalent to the 
commonly used R — K>5, R — K> 5.3, and I ~ K > A. 
In the survey, reaching approximately K = 17.5, we find 
54 EROs. The area covered is 2850 arcmin^. 

Taking advantage of overlapping mid-IR data, we 
search for dusty EROs, since only these should be de- 
tected with the used 15/L(m ISOCAM band. Only one 
is found from our conservatively constructed catalogue. 
Taking into account detection limits we limit the number 
of very strong starbursts {SFR > 2OOM0year~^) in the 
bright K <17 - 17.5 ERO population to < 10%. 

We also make use of a J—K vs. optical-infrared colour- 
colour diagramme to separate EROs, and find that the 
fraction of dusty ERO population is < 10 — 40%, depend- 
ing on the colour used. There are more dusty galaxies in 
the i — K based ERO selection than if r — if is used. 
HYPERZ photometric redshifts and template fits are also 
employed: nearly all redshifts are in the range z = 0.6 — 1.8 
with a strong peak at z ~ 0.8. Approximately 90% of the 
best-fit SEDs are those of evolved stellar populations. 

We find a considerable amount, ~ 60% of all our EROs, 

r — K > 5.5 EROs which are not EROs with the 

1 ~ K > 4.4 criterion. Using models we interpret these to 
be early type galaxies at redshift of z 0.7—1.1. They are 
interpreted to be the counterparts of local 2 — 3L* galax- 
ies, and their resulting space density is approximately 

2 X 10-5 Mpc"^ 

Cumulative number counts are provided for the EROs, 
extending the available ERO counts to brighter limits 
than previously. Our counts are consistent with litera- 
ture counts in the overlapping magnitude region with same 
colour cut-offs. 

Our ERO number counts, as well as other literature 
data, are well fit by pure luminosity evolution models. 
Formation redshifts for early type galaxies in excess of 
z = 2.5 are required to fit the ERO counts, and z w 3 
is favoured. However, the range in the colours of EROs 
suggests also a wide range in formation redshifts. 
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Table 1. Photometry of the whole ERO sample. Sources 1-17 are those with r — K > 5.5 and i — K > 4.4, while 
18-50 and 51-54 are selected by only r — K > 5.5 and i — K < 4.4, respectively. 



14 



Vaisanen & Johansson: Bright EROs and massive galaxies at z ~ 1 



Sample 


N 
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<i~ K > 
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(9) (10) 


(11) (12) 
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6618 
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17.5 


3.91 


3.09 


50 


21 17 


0.007-0.018 0.004-0.007 



Table 2. Columns (2) to (5) give the total numbers of sources and stellar objects, the surveyed area, and the average 
limiting magnitude. The rest of the columns refer to galaxies only: columns (6) to (7) give mean colours, though note 
that these are lower limits since optical non-detections we ignored. Column (8) gives the number of EROs selected by 
r — K > 5.5, (9) by j — if > 4.4, and (10) by their combination; column (11) is the range of surface densities resulting 
from either r or i band based selection; column (12) is the range of cumulative surface densities at if < 17 
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Table 3. The sample of Mt. Hopkins EROs. Cumulative counts are given for different ERO colour selection criteria. 
Total numbers (N), Fraction of EROs compared to the total galaxy count ("Frac") is calculated using galaxy counts 
in these same fields (Vaisanen et al l2()()()|l . At if < 17.5 N shows the completeness corrected value, and the raw count 
is given in parentheses. 



